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Abstract of JP2002201998 

PROBLEM TO BE SOLVED: To improve the 
accuracy in operating the amount of the air 
charged in a cylinder in transition. SOLUTION: 
An opening command value is determined 
corresponding to an accelarator operation 
amount or the like during the operation of an 
engine, and the output timing of the opening 
command value is delayed by a predetermined 
delay time. An estimated variation of a throttle 
opening is operated by an electronic throttle 
valve on the basis of the opening command value 
before delay, the estimated variation is added to 
the present throttle opening (the output of the 
throttle opening sensor) to determine an 
estimated throttle opening at an intake valve 
closing timing. Then the estimated amount of the 
air charged in the cylinder is temporarily operated 
by the intake system model on the basis of the 
estimated throttle opening, and the obtained 
amount is differentiated and integrated to operate 
the estimated variation of the amount of the air 
charged in the cylinder. The estimated variation 
is added to the amount of the air charged in a 
base cylinder operated by a base intake system 
model to determine the final estimated amount of 
the air charged in the cylinder (the amount of the 
air in the cylinder determined at the intake valve 
closing timing). 
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/Pa *^5y-ji4-r*"f— •^Mf^mmT^tA 

f (Pm/Pa ) <Ot—7MZ. Wckm*l £LT 
jEMitTZt. 06K^TJ:^)5:ft«r*sn-5„ C© 
f (Pm /Pa ) «. JJB (3) 5£ s (4)3&>€>9j& 

*»fcJ:9tt:. &ffltotett:ft©ffiK:&&&i,>;te«>. 06 cd 

Wt?tt. Pm/Pa> 1©4#«C. f (Pm/Pa) = 

[0047] Pm /Pa > l©i#tcf (Pm 

/Pa) = 0K"rS£. Pm / ¥ <dfi\ W&V&mt % 

IttUa^ h^iXSMIGin. ^ffi'MMEEPm . ^ 
8'J)gfrt^«^fiGcf) jWSIo L"Cr»>* 
Sttrtl***. COJHHtt. Pn /Pa#UflSE©«« 
"Cf (Pm/Pa) ©SEffc****<<C*Cii. ftft 
I^HERltcPn /PajWRflLh 1 filhCC&S&c f (Pm 
/Pa) ftSO-Ctf-KStl&fcti). ]©mwjlfeB$© f 
(Pm /Pa ) ©^ffc*CF*«lJtCfcSfc#t?**. 
[0048] CW*mt Lt, ( 1 ) Ttt, 

f (Pm /Pa ) ©f-^^iaTK^TJrOtClS^L/ 



9< Pm/Pa<l©£#tCf (Pm /Pa ) =IE 

30 ©<B 

Pm /Pa = 1 <D±%tCf (Pm/Pa ) =0 
Pm /Pa >l©i^tCf (Pm/Pa ) = fi©fil 
CCR£t/?t>«. C*ifc«fc9. f (P«i/Pa)Of- 
^JH*. Pm /Pa = 1 **KbT±3MKlK-j-4*tffc& 

[0 049] H7©<fc5£^fc#t4©f. (Pm/Pa ) 
Of-^*fl!l>4 4. Pm/Pa*5 1«ja-C^8tir-5. 
iWHwilfefStc. f (Pm /Pa ) ©gf fcflWMWfc <> 
©£ft5„ C©ft:«>. IftfUfi-^^oaiffffi (*P* h 
40 JUffia^mGinXB^S'J&JftEEPm Xtt^TOWrt^* 
SSlSGcf) tfJWtriCir. H9fCwTJ:5»c/ 

jSGcf) *^<fc3-a-sc^*!-cir. ^>^>* r *esrjh 
[0050] ^ao* n » h ;i/iiaasg^,a?amsi5K: A^j-r 

«BURBEPmC1-i3W]il,»6*i**». ®fUHz>1H 6© 
[0 0 5 1 ] i?c. ^S'jXU-v h/biMiS^tfiGin©?* 



(10) 



17 
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mm&m ( i ) t». immMZffimitrzfc&ic. m 

[0 05 2] JVC. TOJftfitEPm &*Wttl*3aSK£« 

l l ©*»PSt?©ftfR»K (fiTF rxa 7 
KS«i?©ffi|ij*®M-rs<i:. ;x© (5) sSrgtsns *io 

Gcf=7? • (Ne /2 ) • Vc 

n : 

Ne : x>y>01gaa! 

Vc : 

v IM : *n * v >vTmmms&<Dttgm 

[0 054] ccr, {*»8rt57? tt. mA&sm&icj: 
^S'J®ME©Ull©fiiPmCi-l)-C*.5> <> 20 



d/d t • Qm =Gin-Gcf 



(5) 



CCT, Qm v h;uTS£®maif&F , 3©^Jt d 
/dt • Qm IWO'j h;L-T*eR»aSSF^©^MS©^ 
ftS. Glni*F»*u» h*9&£4Ut GdttflBlli 

[0 053] JbT^tR»afStcmft©t^ 

@*fis*affl-rsi. ^© (6) ^fca-rw«*5»6 

ft£„ 

• (Qm /V IB ) (6) 

* [0 055] «^*^1/©*^b#j8Rt ib h: 

JFSD (7) St?«§n5.- 

T I( , = 2 • V IB / (Vc • 7j • Ne ) (7) 

±ia (5) ~ (7) 3&6>&#c© (8) s^^m^n 

■So 

d/d t • Qm =Gin-Qm /x lu (8) 

[ o o 5 6 ] ±gg ( 8 ) s«ag«5sreftsfc». cm 

«-?*U®)J.-y h 2 5Tig#&SlT-£5J:Mc^£;fci* 



77= f (Ne . Pen ) *• 

{Qm(i)-Qm(i-1)} /Ts = Gin(i) - Qm(i-l)/r I1( (9) 

ccr. Ts >i*9ni-c*&. *tt§tauiigrt©£atSQin ©»jiw*J^©j:5(c**a 

C0 0 5 7] C©(9)S58It5i, XUjW* Stl*. 

Qm(i)= (Gin(i) - Qm(i-l)/r IM } • Ts +Qm(i-1) [kg] 

(10) 

[0 05 8] Sfc. XD? h*T8lfflWiaiS{cjllft©tt #mfiQm frh^fflmfS&Pm *M^£5£#&©<fc5(C 

Pm =Qm • R • T/V IM [Pa] (11) 

R : ^,<*5£ffc 30^ [ 0 0 5 9 ] ±13 ( 1 1 ) j£t ( 6 ) 3<b*>&. #© 

(12) srasna^wiwwai^WGcftoiSMac 



f«UR*?*©^!RjnEEilW*ti, ±12 ( 1 0 ) 5£<b 
(11) SSffl^t, ^JGR^EEPm ^iHSTS, ♦ 
Gcf= 77 • Vc • Pm / ( 2 • R 

«a^*^©^wiBfrt*aiffiaw[jitspB:. ±ta < 1 
2) sc*m>r. {g©^»jiSiF«95t«^msGcf*?iiff-r 
s. 

[ 0 0 6 0 ] m 2 icm? <fc 5 &t\ ©SUf^e-r;!/©^;*? 
(^©^S'JtgirtJtfil^MSGcf) fit. Sfc#g#t (d/d 

t) fcA^sn. ^^'^mBt s ng©:g##*«> 
6n. *©M3-#«a«3i ( J ) -eat^sft*. -e©n 

fl^fSW:, M*«*fST A U©Stl? # Y 5 > ^ <tSrt=fc 
tt^eig©^*^ 5>^) friir&MJW7ffl2J$> 
^^■C©Hf^Tin.i "CiS. «#SS?t ( S ) -C»#Lfc 

fiiis . r*/ t-k # w = > ^ * -e©®rt?t«^*s© 

AGc *^-*!8«»*f , *fcJ:g»J|£LA:^-;*tti*j 
3fe«S«fiGbase{C»PfiLLr. S^&^S'»3fe»S 

SIGc (^^:/&i*^s>^r^ir3t£rt^ 



T) [kgr/r e v] (1 2) 

[0 06 1 ] #tc\ -«-^«|*K(E«ffi««©S!»^«* 

# 1 4©m^j vRAmmmm) <<cm-?^xmnhtdm± 

©fgff*3^t*^fir*-5„ fie^r. c©^-.xfSf*#ti* 

rfctcj;sfSrt^«^««©^b»*^*n^o. -m 

{c % i77ny-$ l 4©&;ft*>6fifF f 33SW£m£3:?* 

tmm&ifi&&ifi. a?gB#ic». xT7n>-5fi4 

©JSSiin (Mil«Oi77D>-$ l 4©1^ 
tt, ir^ay-* 1 4©-fe>-^gPi#©b- hv^tc 
iS^ilti) *s#?4TSfc». jffiam©J£9(feMBC» 

so [0 062] fcr. *HSS0fig ( 1 ) T«, il?gB#© 



(11) 

•^-Xg}rt3S*^mSGbase?r?||ff-r-S.. CO-^-^f 

Gbase= 1 / ( 1 + T,„ • s ) • Gdlay 
Gbase: ^-XfUrt^il^Mffi 
Gdlay: lD§jiniffitg3l©ffi;*J 
r„ : B$3£$fc 10 
[006 3] CO^-^Km^^rt'^i^ifcT^tt, 

r,. = 2 • V tM / (Vc • 7? • Ne ) 

Vc : ZsV^fgffi 

v : imwm 

Ne : x>y>0^3IS 

x>i^>[5lg3iaNe iKMSP (®MEE-fe>lM 6 © 20 

[0 06 4] CCD^^i&^-^tR^^&^KCj;^?^ 
Ufc^-^1giP*g5t*SmmGbasei, ^SflXP? h*H 
^*>6^rL/cfgfF^*£MS©f^»bfiAGc £ 

srafit^r. ^t^ft^iuitn^ix^sGc (crsia 

[006 5] «±Ui?gUfcBI2 (D#^d ? *©&fjt«\ 
01 075MS1 8©&JU-^>fc:J:-3 Testis. 30 

[0 06 6] [y-f SI 0<DM-<(Z/)\s-9- 

oot, ®m-rzm i layxu-j b)ifflmm>i-*> 

total =&3f^HTdly jiM3i±&XP v h;I/S5i$0^ 

£3£frr&. Co®, o otcit*'. ?&i£-r& 

[oo67] co^ x^?^ o otcji^. 

»3«2«#fSTp C©&. ^f-^40 0 

U fi#^lfrfc#^S*K£MjB&&Kload (ftliSiilfiE 

mo&mvmiEm&K c*s*pgs=fsT P team o so 
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[0 06 8] [xp h;l/jffl^f(|^l'-5 L >] S 1 1 © 
x p ? h ;i/il5i$iJfflJ^-^ >«> 0 1 o <dj A 
><D^f o or-^fsn^^^;i'-^>-r-*5„ 

r. resJHMi* total fctssrr.*.,, c©bl wsm^ 

fit* total ». T*teJ\smft&l l CJfcttcm?fcffl&<t>peda 

u, 7-f FJHsWEiss/gu® (isc) tccts^gas 

<t> total = <t> pedal + isc 

[0 06 9] C<DX^-;-7 1 0 1 ©MS^ftlt*©® 
[0 07 0] C©lt. Xf-y-fl 0 2fCit*> XPy h 

ate**, xu, HJi/jMaww»±*ff«cFiaat4ft*. 

[007 1 ] ^Kii^^KiKft}*. **WK:x>y>@ 

/c, 7-/ kjwwei&b. r-f FJna^iiaMai ( i s 

C) *ifPlftL-CT^ F*|g(gatflEt7 - K/*» 

aS*iJ@l3^XP 5. F;HSM»if^tt7-f F^lllg 

-h^;H^©i^B. ai£#^tx— »if (x>^> 
©£P* D 5 DJtglSJiS* -5 tcti>. - * - t- 5 Jl/tfclS 
©i#tc. XPy h;l«SW»*tf5i. U—>>^* 
^B${Ci>y>lie2S©lLi:*s>lAianr. ate^fc 

[0072] -g-C-C, 2(c||S60SI ( 1 ) -C«> XP 9 V 

p » h *^ifc-r ici-c, x p y h »«tl 
[0 07 3] feU, Xf^l 0 2-C. XP-y hJUiljg 

<«•?) ©Wfldfi^ffi* total Ci)*iBlSl±-r«:*-* 
[0 074]— ^r. ?yi02T. XP^H;It1JI 

JW^©«tS«: <£ o t . -X© J: ^ or x p ^ h ^iSSfiJ 



(12) 

21 

total (DilfilSHTdly £j&rr£„ C©l£. 3JP$R3 

SlXA-t^W * ^ 3 > -COBSIffl T i m #> hna-?. P 9 
h;l^*^A©M*B#ffflTth£gL?|l>fcB#W (Tdly 
= Tini -Tth) tcSCStiS. flu *P4«S«TA 
ucomn 2 4 3 > y*» fcfWfW * 3 > ^2: -c© 
BSfflTirvj #*ttW*BT-thJ:D fcg< (Tin.i 
-Tth<0©f§i^) B. «BPSIBTdly*0 
[0 0 7 5] C©f£, Xf^l05 tcjlt*. jB&$m 10 
Tdly |*3©-?>-7-</ >^Cdly *&£(C «fc 0 MtW 

Cdly =Tdly /Ts 

cct, Ts ';>^p H i-c&.2>., 

[0 07 6] C©f£, Xfj^lOB fcit<&. 9 
feil^fflTdiy rtcDif^^y >^SiCdiy micmM-U 
tcmS.tB^U<t> total (i-Cdly) £|gSblHlS§ 3 2 (C 

Httrr*. c*uc<fc9> ^SJi^iii* total ©ffl2>*-f 
5>**»WIBTdly MttaBBd-ttS. Cti^f? 
^"1 0 2~1 0 BCDteS&S^F^Oififfl-CC^f 1 ^ 1/ 20 

[0 07 7] [ ^HSfF^ift^MfiSI?^ - ] H 1 

2©?a<jjiji^«s^«*M'--^>B. si i©y-r 

>Jl< - >©*■?- 9 :/2 0 OfUff^n-S.lf -7>\,-*l/ 

[0 0 7 8] ^-^asjiasbsnsi. ifxf^ 

201T, 1 3©^S!fffi^ffi||#^-^>?: 

^jWf^UEEPm (mM'i>l>7mZJ 5>i«Oft 
fiflE) *?SISfrSo C©f&. 0 2(Cii* > ^ 30 

W««]£P« iffllit. gSttcJcQTMRftaEK'SflUI 
Gcf(i) feffiff-TS. 

Gcf(i) =77 • Vc • Pm / (2 • R • T) [kg/ 

rev] 

77 : ttflltt* 

Vc : S^y >#S» 

r : imim. 

T : «R»SS 

[007 9 ] C©f£. ^2 0 3 Kjt.^ El 1 1 © 
^^102 £|5l C^ffi-c. * n 9 h ;vaSM^ih 40 

jtm»WS©«S[^ 5 > feB&CWT'm Z J 3 > 
y $ T©«roisa^Mt(O^SBEifc«A Gc £:0 (CfS5£ 

[0 08 0] xn, h;uS52$U©^±^*s^BS 
£©*§£«, 0 1 \<oxuv h)\>W£MM)\>- : ?'Zslc£. 
ot^n 9 h;bjlJi^€:l©Sbr, Xf? ^2 0 3* 
6X^5.^20 5^cil*. ^S'JB^MTini F*3©1f>^'; 50 
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Cp = Tin.i /Ts 

CC-C, ^SiJ^WTini B. IK**»gS«©SiJ?* >C 3 > 
^p.Ka^T'H^'f a>y*-c?©ii*ffl"C**>. Ts 
B-9->7">J>^IWcAS 0 

[0 08 1 ] C©f&. XT" 57*206 {Cit*. M£W*f 

i©ig * >f 5 > y*> e>©«A- ji, 7^ -/ * > y s t o 

fSF*33t«^SGcf©^I>J^taAGc ^^KcfcOjg 
AGc = Gcf(i) -Gcf(i-Cp) 

cc-e. Gcfci) b, ^<o*mmp&e*smM. (-ot 
k> wi&iZAxis if(omn % -r s > y©^a(iMrt 
-c* o . Gcf(i-c P )B. m&£ <o h^mm 
faTini >y«cp MicmnLtcWipffim 

{^<o tmmmmvmm- * ■< 3 > ^©ns^a 

[0 082] ^S^fbSA Gc ©iftHf£. X f 9 7*2 0 
7fcii«&. '<>-XWW5t«^MM?llfit-'l'"5=-> (lit 
T) *|gff U ^-X(gfrt^«3gM«Gbase^jS»t- 
-S„ C©l^> I7?ay-* 1 4©ffl^©JDE^iin^I& 

Kite j: o ^-^isrt^a^amGbasG^vgm-rs, 

Gbase= 1/(1+ t la * S ) • Gdlay 

^/cAfc. ^-xfgsrt5t«sma©M@s^iii^^-ca 

L/-Tl^*s, m^JtSlJ.- •> h 2 5B. ±^i»tS{{tL 
7*cS?:fflt^-C^-^Si^«^SGbase5:j||#-rS. 
[0083]-eoi. 9 -^2 0 8^11*, t©^- 
^K^JK^fiGbaseK:^ ±IB^7- ^2 06 -C** 

MMGc *&tbZ>. 
Gc =Gbase+AGc 

[0 084] [?®MSg,1EffiW>l<-*>l El 1 3 ©^«>1 
6R^E«gtJl/-^>B. IS 1 2 ©^aiJ©rt^«^«a?S5 
»Jl/-^>©^-r-9^"2 0 1 f Hff 3nS1J-^I/-5 : -> 
^-^V^iajidstiSi. *T^f 9^2 1 

l r . §114 ©^jfflJ^ n 9 h 

gE-r^o c©f^, xt-^2 1 2icjt^. ^-r^si 

6©®«^*7 r ^S^ff>'l'-^>**tfL/. Rmifi 
^^©^e^i/^^rL^MSL-rs. c©^, ^-^9 
7*2 1 3(Cjl*. *d9 h ^TSE©mjlSgrt©^mMQ 

Qm(i)= {Gin(i) - QmCi-l)/T,« } -Ts+QmCi- 
1) 

CC-C. Qm(i)B4-0©* P 9 h ^T^K^a^F*?©^ 

ma. Qn»ci-i)BWing©*n9 F-^T^Kmasgrt©^ 

[0 08 5 ] C©*£. Xf-9^2 1 4^*, 7LU9 h 
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Pm = Qm • R • T/V,, 
CCt. RS*fW££§&. TttBWtfiflE. 

[0 08 6] C©f&, Xf9^2 1 5tCit#. ^-lelCD^- 
ffiflftflS P m(i ) i SfrlUCD^aiJK^IEE P m(i -1) i (D¥-i$)m 

Pm CO = {Pm(i)+ PmO'-l)} /2 

[0 08 7] [ toj* p ? h ;UffljS^ma?H^L;u- ? io 

> ] H 1 4 <DfWX a v h JHMjftS«S»ff * > * 

Pa 



Gin = #* • A • 

Gin : 
U 
A 
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n-st. mf»^2 2ir, §aw*H.i 5©^»j 
A s >i/©fffl^a y h;ngge f £itg-rs. c© 

fit Xr-^7'2 2 2&Cjt<*K *^EEPa. S&gUfiKTK. 
Olg EJ©f-«'M*JE P mCi-l)fc^#jAif . 
[0 08 8] £©*§. -?7'2 2 3 JCJl*, 

[0 08 9] 
[&4] 
f(Pm/Pa) 



^nyhywaiS^SUI[kg/sec] 



Pa : zfc&JECPa] 

Pm : ffi^tflEtPa] 
R : 

T : K^CS^CK] 

fCPm/Pa) : 85SLEPml:^:^(tBEPa«3lfc-edt*6«lffiffi 



[0090] C©^, M • Ait. "miXUy hJl/Hge 
f ^^-^itST-^^ffil/. f (Pm/ 
Pa ) B, Pm /Pa # £T5EI 7 ©?— ^ 

A*6J|Effi-rS. BMEPmiJ. wTlil©^SiJ5R»BEPm 30 
(i-l)>CP*fflt,»6n. AMBEPai^^aSTK. -g-ft«eft 

■b>v<D^tamami>^ti^. ft, ^MEEPa(*amA 

[0 09 1 ] [ffxn-; hJVHflE«»l'-^>] M 1 
5 © ^SJX n 9 Y Jl-PBflSBflUl'- HI 1 4 ©^S'J 

XD^l- -U/jfit^afiSS-'l'- >©^f- 

^tf 3n5-9-^-^>f*»), ^pit«*©siH-ei^ 

h^Sja?S@iLT©!9«lI?rm/c-r 6 
[o 09 2 ] ^v-^^iasd^n^i. If^f^ 

2 3 17?. T^-fe^JSflFfi^Cf&DT^aF^ffi^tota 40 
1 *S5c-r?)o C©KS. gl^g^fii* total 
^&WS5£)£D/cH#|gS<*>peclal «t. T-f F;HH3fciI 
gftj® (ISC) tC<fc£g*Mg<*>isc ^©JH^&^jfS 

</> total = <p pedal + <t> i sc 

[0 09 3] COfi. 7^2 3 2 Kit*, XO-^ 
J^S-fe>^l 8-C^i±lLfc^aE©Xay hfrm&dZ: 
U*tiLA,tc&. Xf 1 » 7*2 3 3 Kit*, 04 KTjVTS-?- 

fHamtt£tc£o?^tt0M^4tf* total so 



(f§fP^ti*^mfi©^I'J# W 5 >^> *»6«flV<Jl'^H 
3>^^-C©B#raTTn.i Y^mSL<0 : fmSt 

ftat-T**. eu (Rsuw^ih*^ 5>^s-r?©i$fffl 

Tini b;i/^^f-A©^|S;B#raTthJ:0 4> 

C©iSE|*^raTthrt©xay h^UBS©^ 

^m<t&A e ft £ . 

[0 0 9 4] C©M. Xf»^2 3 4tCl^ 3HSE©* 

^ (xtt*w«raTtiisa«) ©^aj^n? h^Hijjtr 
[0095] [ msB&^ftotimm&m-fr-* > ] h 1 

e©«WUR- ! ef s Jl*3eSB!{ML>l'*-^>tt. S13©^IiJ 
!RaESfjWf>©^f^2 1 2-C»f$fti-^y 

v?24i aas-r*H 1 7 ©^^m#^-^> 

*SttTLr. #8B»*r?ft«(W-rS. C©». XT'-;? 
2 4 2 tCjt^ v «fJUI»aR t : « Sr^K «t «5 ^gC-T 
4. 

r 1 . = 2 • Vx./ (V c -Tj-Ne/60) 
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ffi) . V c (*->'J>^^a . Ne tfri>i?> 

HIIEaS ( r pm) -C£>£„ 

[0 09 6] [»a«)3S8»JU-^>] HI 7©#«» 

*;U-?>^l6?n-5<!;, 5 1 

■C, fu0©©S*JEPm(-i-l). AMEPa. (KMSJKT. 

^cSTHW^r^jitf,, C ©f£, Xt- y 7 1 5 2 Kit 10 
^. Pm/Pa. x>i>>[IieitgNe . 5 

[0 09 7] -[W«aiiE^-^>] HI 8 ©«***§ 
iEJl'-5 1 >«. H 1 0<DJ JZsJls—T-XDZ.?-? 7*4 0 

or»f3n-5-9-^-^>-c*o. m&nuamn' 
WfiM»#Si or©!asw*»fc-r. 20 

[0 09 8] *il/-^>*i*ag(|$nSi. S^T^ 
40 It?, r^-feA»ffKJ:*-ft-fl80ft 

0 2 jtns&to (is^a^fio^tt) {c^-r 

S^£fME^Kload£/h3tHSKi fcSST*. C© 

se&». *£i«gtt ( 1 ) (Dtsftimzm&cDmnysm 

[0 09 9] T^-fe^f^fciSftjw^Sjj-c^li 

^:E{fc&Kload£;ft:i*H>filK2 {Cta^T.?,. COlft 

[0 10 0] fedOi^Kl/t, ^-^4 02XtJ4 
0 3T. fi-^S&KftTSj^ffilE^&Kloadfc^l, 

v74o 45CJ1*. fi^atit'^©saK^ 
Ts&a©*gm«E«£f*Kc (m«£**it7 -f - f-k 

7^MCE0MK. ^sntiEflttt«) 

fi^Kioad. kc R.vmmmmmmTv *iur, 

Wttt»Wtt« (*^>VWXl»§> T A U*;*5»c J: «3 fH 50 
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TAU = Tp xKloadXKc + Tv 
[0101] JM±»qi^-t«*--f>KJ:r>raiH/fc 

£H 1 9©£-f + - FtC^bT^S. i>^>®^ 
rf -fe^»ff«*K:ifccriBa3B^ffl* total * 

lags u . c (om&iu^m <*> total ©m* z a s > 

SB-JHTdlv 7ttt«ES-tfS. COR. agBfflBTdly 
tJ> H3{C7K-r<i: i 5(c. j^Wg^STAU©^^ S 

>^ (f»ft?c«£«s©^«'j*-r s>so a>&<R£vw 

5>^*-C©^Tini Ti^STXP ? h^^> 
XT-A©^i*B#raTth^L"5K»fcB$fa (Tdly =Tin 
1 -Tth) censes {BU jMS*fSTAU©?|| 

irvj SfXO» V)lisX^J*<Dm$tP$fflTthJ:K)i> 
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(57)Abstract: 

PROBLEM TO BE SOLVED: To improve the accuracy in 
operating the amount of the air charged in a cylinder in 
transition. 

SOLUTION: An opening command value is determined 
corresponding to an accelarator operation amount or the like 
during the operation of an engine, and the output timing of the 
opening command value is delayed by a predetermined delay 
time. An estimated variation of a throttle opening is operated by 
an electronic throttle valve on the basis of the opening 
command value before delay, the estimated variation is added 
to the present throttle opening (the output of the throttle 
opening sensor) to determine an estimated throttle opening at 
an intake valve closing timing. Then the estimated amount of 
the air charged in the cylinder is temporarily operated by the 
intake system model on the basis of the estimated throttle 
opening, and the obtained amount is differentiated and 
integrated to operate the estimated variation of the amount of 
the air charged in the cylinder. The estimated variation is added 
to the amount of the air charged in a base cylinder operated by 
a base intake system model to determine the final estimated amount of the air charged 
(the amount of the air in the cylinder determined at the intake valve closing timing). 
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CLAIMS 



[Claim(s)] 

[Claim 1] In the internal combustion engine having the electronic throttle system which drives a throttle 
valve with a throttle actuator and controls a throttle opening An opening command value operation means to 
calculate an opening command value based on an accelerator control input etc., A delay means to delay the 
timing which outputs the opening command value calculated with said opening command value operation 
means to said throttle actuator, A throttle opening prediction means to predict a subsequent throttle opening 
before the delay output of this opening command value based on the opening command value before making 
it delayed with said delay means, and said electronic throttle system response delay property, A restoration 
air content prediction means in a cylinder to predict the restoration air content in a cylinder based on the 
throttle opening predicted with said throttle opening prediction means, The control unit of the internal 
combustion engine characterized by having a fuel-oil-consumption operation means to calculate fuel oil 
consumption based on the restoration air content in a cylinder predicted with said restoration air content 
prediction means in a cylinder. 

[Claim 2] Said restoration air content prediction means in a cylinder is the control unit of the internal 
combustion engine according to claim 1 characterized by predicting the variation of the restoration air 
content in a cylinder to intake valve close timing based on the throttle opening predicted with said throttle 
opening prediction means, adding to the restoration air content in a base cylinder which calculated this 
variation based on the current operation parameter, and predicting the restoration air content in a cylinder. 
[Claim 3] Said restoration air content prediction means in a cylinder is the control unit of the internal 
combustion engine according to claim 2 characterized by to predict the variation of the restoration air 
content in a cylinder to intake valve close timing using the inhalation-of-air system model which applied the 
law of mass conservation to the inhalation air which considers that throttle opening which inhalation air 
passes is an orifice, and flows a throttle passage air content and a throttle down-stream path in integrating 
the variation of the output of this inhalation-of-air system model to intake valve close timing. 
[Claim 4] The formula which calculates a throttle passage air content with said inhalation-of-air system 
model is [Equation 1], 
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In case it is alike, and is set up and said restoration air content prediction means in a cylinder calculates said 
throttle passage air content, f (Pm/Pa) is Pm/Pa. It is the control unit of the internal combustion engine 
according to claim 3 characterized by computing from the table made into a parameter and computing mu-A 
from the table which makes a throttle opening a parameter. 

[Claim 5] The table of f (Pm/Pa) used for the operation of said throttle passage air content It is set as an f 
(Pm/Pa) = negative value at the time of f(Pm/Pa) =0Pm/Pa >1 at the time of f(Pm/Pa) = forward value 
Pm/Pa =1 at the time of Pm/Pa <1. Said restoration air content prediction means in a cylinder is the control 
unit of the internal combustion engine according to claim 4 characterized by having a means to equalize the 
operation value of said inhalation-of-air system model. 

[Claim 6] Said delay means is the control unit of the internal combustion engine according to claim 1 to 5 
characterized by setting the time delay of an opening command value as the time amount which deducted 
the dead time of said electronic throttle system from the time amount from the operation timing of the fuel 
oil consumption of a certain cylinder to the intake valve close timing of the cylinder concerned. 
[Claim 7] Said delay means is the control unit of the internal combustion engine according to claim 1 to 6 
characterized by outputting without delaying said opening command value when the time amount from the 
operation timing of the fuel oil consumption of a certain cylinder to the intake valve close timing of the 
cylinder concerned becomes shorter than the dead time of said electronic throttle system. 
[Claim 8] Said delay means is the control unit of the internal combustion engine according to claim 1 to 7 
characterized by outputting without delaying said opening command value when an automatic transmission 
corresponds to either of the neutral conditions at the inside of the predetermined time immediately after start 
up, and the time of idle operation at the time of start up. 

[Claim 9] Said throttle opening prediction means is the control unit of the internal combustion engine 
according to claim 1 to 8 characterized by predicting the throttle opening after an opening command value 
delay output using the electronic throttle model containing the lag unit and rate limiter more than primary 
[ which considers the opening command value before making it delayed with said delay means as an input ]. 

[Claim 10] Said throttle opening prediction means is the control unit of the internal combustion engine 
according to claim 1 to 9 characterized by predicting the variation of the throttle opening to intake valve 
close timing using said electronic throttle model, adding this variation to a current throttle opening, and 
predicting the throttle opening of intake valve close timing. 

[Claim 11] Said fuel-oil-consumption operation means is the control unit of the internal combustion engine 
according to claim 1 to 10 characterized by having a means to amend fuel oil consumption according to 
operational status, and switching the correction factor to fuel oil consumption by the time of the load effect 
by accelerator actuation, and the case of being other. 

[Claim 12] In the internal combustion engine having the electronic throttle system which drives a throttle 
valve with a throttle actuator and controls a throttle opening An opening command value operation means to 
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calculate an opening command value based on an accelerator control input etc., A throttle opening 
prediction means to predict the throttle opening of intake valve close timing based on the opening command 
value calculated with said opening command value operation means, and said electronic throttle system 
response delay property, A restoration air content prediction means in a cylinder to predict the restoration air 
content in a cylinder based on the throttle opening predicted with said throttle opening prediction means, 
The control unit of the internal combustion engine characterized by having a fuel-oil-consumption operation 
means to calculate fuel oil consumption based on the restoration air content in a cylinder predicted with said 
restoration air content prediction means in a cylinder. 

[Claim 13] Said restoration air content prediction means in a cylinder is the control unit of the internal 
combustion engine according to claim 12 characterized by predicting the variation of the restoration air 
content in a cylinder to intake valve close timing based on the throttle opening predicted with said throttle 
opening prediction means, adding to the restoration air content in a base cylinder which calculated this 
variation based on the current operation parameter, and predicting the restoration air content in a cylinder. 
[Claim 14] A restoration air content operation means in a base cylinder to calculate the restoration air 
content in a base cylinder based on a current operation parameter, The inhalation-of-air system model which 
applied the law of mass conservation to the inhalation air which considers that throttle opening which 
inhalation air passes is an orifice, and flows a throttle passage air content and a throttle down-stream path is 
used. A variation prediction means to predict the variation of the restoration air content in a cylinder to 
intake valve close timing based on the variation of the output of this inhalation-of-air system model, A 
restoration air content prediction means in a cylinder to add the variation predicted with said variation 
prediction means to the restoration air content in a base cylinder calculated with said restoration air content 
operation means in a base cylinder, and to predict the restoration air content in a cylinder, The control unit 
of the internal combustion engine characterized by having a fuel-oil-consumption operation means to 
calculate fuel oil consumption based on the restoration air content in a cylinder predicted with said 
restoration air content prediction means in a cylinder. 

[Claim 15] An intake-air- flow detection means to detect the flow rate of the inhalation air which flows an 
internal combustion engine's inhalation-of-air path, The inhalation-of-air system model which simulated the 
behavior of inhalation air until the inhalation air which passed the throttle valve flows in a cylinder is used. 
It has an operation means to input the output of said intake-air-flow detection means into this inhalation-of- 
air system model, and to calculate the restoration air content in a cylinder which is the output of this 
inhalation-of-air system model. The time constant of said inhalation-of-air system model The control unit of 
the internal combustion engine characterized by setting air content change as a small value so that a actual 
twist may also appear at an early stage. 

[Claim 1 6] A means to presume the current restoration air content in a cylinder based on a current throttle 
opening, A throttle opening prediction means to predict a future throttle opening, and a means to predict the 
fUture restoration air content in a cylinder based on said future throttle opening, A means to add to the 
restoration air content in a base cylinder which calculated the deflection of said future restoration air content 
in a cylinder, and said current restoration air content in a cylinder based on the current operation parameter, 
and to calculate the final restoration air content in a prediction cylinder, The control unit of the internal 
combustion engine characterized by having a fuel-oil-consumption operation means to calculate fuel oil 
consumption based on said final restoration air content in a prediction cylinder. 



[Translation done.] 
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* NOTICES * 

JPO and NCI PI are not responsible for any 
damages caused by the use of this translation. 

1 .This document has been translated by computer. So the translation may not reflect the original precisely. 
2.**** shows the word which can not be translated. 
3. In the drawings, any words are not translated. 



DETAILED DESCRIPTION 



[Detailed Description of the Invention] 
[0001] 

[Field of the Invention] This invention relates to the control unit of the internal combustion engine which 
improved the operation approach of an internal combustion engine's restoration air content in a cylinder. 
[0002] 

[Description of the Prior Art] In order to correspond to the emission gas clarification regulation which 
becomes severer future still, it is necessary to high-degree-of-accuracy-ize Air Fuel Ratio Control (fuel- 
injection control), to calculate the air content (restoration air content in a cylinder) with which it fills up in 
an engine cylinder with a sufficient precision for that purpose, and to set up the proper fuel oil consumption 
corresponding to the restoration air content in a cylinder. The operation approach of current and the 
restoration air content in a cylinder currently generally used detects an intake air flow with the air flow 
meter installed in the upstream of a throttle valve, and is divided roughly into the method (mass flow 
method) which calculates the restoration air content in a cylinder from the detection value, and the method 
(speed density method) which detects an intake pressure by the intake-pressure sensor, and calculates the 
restoration air content in a cylinder from an intake pressure and an engine speed. 
[0003] 

[Problem(s) to be Solved by the Invention] By the way, the timing which the restoration air content in a 
cylinder decides is intake valve close timing which an intake stroke ends, and the timing which computes 
fuel oil consumption is a front [ timing / intake valve close ] (it is because it is necessary to perform fuel 
injection before intake valve close timing to make an injection fuel inhale in a cylinder). Even if it computes 
the restoration air content in a cylinder by which approach of the conventional mass flow method mentioned 
above and a speed density method, therefore, to a transient The restoration air content in a cylinder changes 
from the calculation timing of fuel oil consumption before intake valve close timing (definite timing of the 
restoration air content in a cylinder). Consequently, the ratio (air-fuel ratio) of the actual restoration air 
content in a cylinder and the fuel quantity which flows in a cylinder shifts from a target air-fuel ratio, and 
there is a fault that the Air Fuel Ratio Control precision of a transient worsens. 

[0004] This invention is made in consideration of such a situation, therefore the object is in offering the 
control unit of the internal combustion engine which can raise the Air Fuel Ratio Control precision of a 
transient. 
[0005] 

[Means for Solving the Problem] In order to attain the above-mentioned object, this invention predicts the 
throttle opening of intake valve close timing (definite timing of the restoration air content in a cylinder), 
predicts the restoration air content in a cylinder based on the prediction throttle opening, and calculates fuel 
oil consumption based on the restoration air content in a prediction cylinder. Here, the change of the 
restoration air content in a cylinder of the reason using a throttle opening as a parameter which predicts the 
restoration air content in a cylinder is mainly because change of a throttle opening is generated as the origin 
and change of the restoration air content in a cylinder can be predicted with a sufficient response from 
change of a throttle opening to a transient. 

[0006] It becomes possible to predict the throttle opening of intake valve close timing (definite timing of the 
restoration air content in a cylinder) by delaying moderately the timing which outputs an opening command 
value to a throttle actuator with a delay means in the internal combustion engine having the electronic 
throttle system which drives a throttle valve with a throttle actuator and controls a throttle opening like 
claim 1 . Under the present circumstances, what is necessary is just to predict a subsequent throttle opening 
before the delay output of this opening command value based on the opening command value before making 
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it delayed with a delay means, and an electronic throttle system response delay property, since response 
delay (dead time) exists in actuation of an electronic throttle system. It can become possible to predict the 
throttle opening of intake valve close timing with a sufficient precision, and to predict the restoration air 
content in a cylinder with a sufficient precision from the prediction throttle opening by this, and the Air Fuel 
Ratio Control precision of a transient can be raised. 

[0007] By the way, the approach of predicting the restoration air content in a cylinder from a prediction 
throttle opening has the advantage that the responsibility of a transient is good, but on the other hand it has 
the inclination for the forecast at the time of a stationary to shift from a actual value according to dispersion 
in an electronic throttle system, aging, a service condition, etc. Moreover, at the time of a stationary, since 
the restoration air content in a cylinder does not change, the restoration air content in a cylinder calculated 
based on current operation parameters (an intake air flow, intake pressure, etc.) is in agreement with the 
restoration air content in a cylinder decided to subsequent intake valve close timing. 

[0008] Then, it is good to predict the variation (a part for transient change) of the restoration air content in a 
cylinder to intake valve close timing based on a prediction throttle opening, to add to the restoration air 
content in a base cylinder which calculated this variation based on the current operation parameter, and to 
predict the restoration air content in a cylinder like claim 2. If it does in this way, it will become possible to 
predict the restoration air content in a cylinder with a sufficient precision by both the time of a stationary, 
and the transient. 

[0009] Furthermore, it is good to predict the variation of the restoration air content in a cylinder to intake 
valve close timing using the inhalation-of-air system model which applied the law of mass conservation to 
the inhalation air which considers like claim 3 that throttle opening which inhalation air passes is an orifice, 
and flows a throttle passage air content and a throttle down-stream path in integrating the variation of the 
output of this inhalation-of-air system model to intake valve close timing. By using such an inhalation-of-air 
system model, the variation of the restoration air content in a cylinder to intake valve close timing can be 
predicted with a sufficient precision by comparatively easy data processing. 

[0010] In this case, the formula which calculates a throttle passage air content with an inhalation-of-air 
system model is good like claim 4 to use a degree type. 
[0011] 
[Equation 2] 

Gin = p m A • ; Pa • f (Pm/Pa) 

Gin : ^n^b/HiiS^SCkg/sec] 

H : assist 

A : ^cr^/VMP^»fB«[[na?] 

Pa : zfc^DPa] 

Pm : «aE[Pa] 
R : 

T : »£«t£[KI 

f (Pm/Pa) : Pm/Pa>l <Dk% f (Pm/Pa) =t^ffl 
A= icr 2 (l-cos2e) 

8 : 2Wi?h/l'VB& 

[0012] In case a throttle passage air content is calculated, f (Pm/Pa) is Pm/Pa. It computes from the table 
made into a parameter, and you may make it compute mu-A from the table which makes a throttle opening a 
parameter. Thereby, data processing of an inhalation-of-air system model type becomes very easy. 
[0013] Like claim 5, furthermore, the table of f (Pm/Pa) It is good to set it as an f(Pm/Pa) = negative value at 
the time of f(Pm/Pa) =0Pm/Pa >1 at the time of f(Pm/Pa) = forward value Pm/Pa =1 at the time of Pm/Pa 
<1, and to equalize the operation value of an inhalation-of-air system model. If f (Pm/Pa) is set to f(Pm/Pa) 
= 0 at the time of Pm/Pa >1 although it does not become a negative value physically so that it may mention 
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later, they are Pm/Pa. There is an inclination which the operation value of an inhalation-of-air system model 
vibrates, and hunting generates at the time of heavy load operation changed in the one neighborhood. This 
cause is Pm/Pa. They are [ that the rate of change of f (Pm/Pa) becomes large in the field of the one 
neighborhood, and ] operation top Pm / Pa at the time of heavy load operation. Since f (Pm/Pa) is guarded 
by 0 whenever it becomes one or more, it is because change of f at the time of heavy load operation (Pm/Pa) 
becomes irregular. 

[0014] As this cure, it is f(Pm/Pa) = negative a value, then Pm/Pa like claim 5 at the time of Pm/Pa >1. At 
the time of heavy load operation changed in the one neighborhood, change of f (Pm/Pa) will become 
regular. For this reason, by equalizing the operation value of an inhalation-of-air system model, the 
operation value of the inhalation-of-air system model at the time of heavy load operation can be stabilized, 
and hunting can be prevented. 

[0015] moreover, time delay Tdly of the opening command value delayed with a delay means like claim 6 
Time amount Tinj from the operation timing (prediction timing of the restoration air content in a cylinder) 
of the fuel oil consumption of a certain cylinder to the intake valve close timing of the cylinder concerned 
from - it is good to set it as the time amount (Tdly =Tinj-Tth) which deducted the dead time Tth of an 
electronic throttle system. It is the time delay Tdly of an opening command value so that the throttle opening 
of intake valve close timing may be made in agreement with a prediction throttle opening, if it does in this 
way. It can set up and the operation of a prediction throttle opening becomes easy. 
[0016] In this case, the dead time Tth of an electronic throttle system is the time amount Tinj from the 
operation timing (prediction timing of the restoration air content in a cylinder) of fuel oil consumption to 
intake valve close timing, although it does not change even if a throttle actuation rate changes. It becomes so 
short that an engine speed becomes a high revolution. For this reason, at the time of a high revolution, it is 
the time amount Tinj from the operation timing of fuel oil consumption to intake valve close timing. There 
is a possibility of becoming shorter than the dead time Tth of an electronic throttle system. 
[0017] This point is taken into consideration and it is the time amount Tinj from the operation timing of fuel 
oil consumption to intake valve close timing like claim 7. When becoming shorter than the dead time Tth of 
an electronic throttle system, it is good to make it output without delaying an opening command value. If it 
does in this way, while not performing useless throttle delay control at the time of a high revolution, the 
throttle responsibility at the time of a high revolution can also be raised. 

[0018] Moreover, like claim 8, when an automatic transmission corresponds to either of the neutral 
conditions in the predetermined time immediately after start up at the time of idle operation at the time of 
start up, it is good to make it output without delaying an opening command value. The time of start up, and 
immediately after start up, essentially, since the engine revolution is unstable, when throttle delay control by 
which an opening command value is delayed is performed, there is a possibility that engine revolution 
fluctuation may become still larger. Moreover, at the time of idle operation, in order for idle rotational-speed 
control (ISC) to operate and to carry out feedback control of the idle rotational speed, when throttle delay 
control is performed, there is a possibility that idle rotational-speed control may interfere with throttle delay 
control, and an idle revolution may become instability. Moreover, since an operator may perform racing (an 
engine carrying out a fore blow) when an automatic transmission is in a neutral condition, when it is in a 
neutral condition and throttle delay control is performed, the start of an engine speed is overdue at the time 
of racing activation, and there is a possibility that accelerator responsibility and acceleration nature may 
impress an operator with it being bad. 

[0019] Therefore, if it is made not to perform throttle delay control like claim 8 when it is the operational 
status in which the adverse effect by throttle delay control appears (at the time [ At the time / At the time of 
start up / of idle operation ] of a neutral), the adverse effect by throttle delay control can be lost. 
[0020] Moreover, in case the throttle opening after the delay output of this opening command value is 
predicted like claim 9 using the opening command value before making it delayed with a delay means, it is 
good to predict a throttle opening using the electronic throttle model containing the lag unit and rate limiter 
more than primary [ which considers the opening command value before making it delayed as an input ]. 
Although it is generally difficult to carry out physical modeling of this at a precision since an electronic 
throttle system has complicated structure By simulating an electronic throttle system response delay 
property with the lag unit more than primary, and simulating the marginal property of the actuation rate of a 
throttle valve by the rate limiter Even if it can build the easy electronic throttle model of data processing and 
does not high-performance-ize especially CPU of an electronic throttle system, the prediction operation of a 
throttle opening is attained. 

[0021] Moreover, a throttle opening forecast may shift from a actual value according to dispersion in an 
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electronic throttle system, aging, a service condition, etc. Then, the variation of the throttle opening to intake 
valve close timing is predicted using an electronic throttle model, and this variation is added to a current 
throttle opening, and you may make it predict the throttle opening of intake valve close timing like claim 10. 
If it does in this way, it will become possible to lessen the throttle opening prediction error by the above- 
mentioned cause, and to predict a throttle opening with a sufficient precision. 

[0022] Moreover, in case fuel oil consumption is amended according to operational status, you may make it 
switch the correction factor to fuel oil consumption like claim 1 1 by the time of the load effect by 
accelerator actuation, and the case of being other. That is, since this invention can predict the restoration air 
content in a cylinder with a sufficient precision to the load effect by accelerator actuation, it can lessen 
amendment to fuel oil consumption. However, since the load effect by the time of shifting an automatic 
transmission from a neutral range to a drive range, power steering, a brake, an air-conditioner, etc. cannot be 
predicted from an accelerator control input, it is desirable to make [ many ] amendment to fuel oil 
consumption to these load effects. 

[0023] Therefore, if the correction factor to fuel oil consumption is switched like claim 1 1 by the time of the 
load effect by accelerator actuation, and the case of being other, amendment of fuel oil consumption can be 
rationalized according to the cause of a load effect. 

[0024] Although invention concerning claims 1-11 explained above was made to perform throttle delay 
control, based on the opening command value and the electronic throttle system response delay property 
which were calculated based on the accelerator control input etc., the throttle opening of intake valve close 
timing is predicted, and you may make it predict the restoration air content in a cylinder like claim 12 based 
on this throttle opening, without performing throttle delay control. Even if such, it can become possible to 
predict a throttle opening using the dead time of an electronic throttle system, and to predict the restoration 
air content in a cylinder with a sufficient precision from the prediction throttle opening, and the Air Fuel 
Ratio Control precision of a transient can be raised. 

[0025] It is good to add to the restoration air content in a base cylinder which predicted the variation (a part 
for transient change) of the restoration air content in a cylinder to intake valve close timing like claim 13 
based on the prediction throttle opening, and calculated this variation also in this case based on the current 
operation parameter, and to predict the restoration air content in a cylinder. If it does in this way, it will 
become possible to predict the restoration air content in a cylinder with a sufficient precision by both the 
time of a stationary, and the transient. 

[0026] Although claims 1-13 explained above apply this invention to an internal combustion engine with an 
electronic throttle system In the case of the mechanical-cable-type throttle system which interlocks a throttle 
opening with accelerator actuation mechanically Like claim 14, while calculating the restoration air content 
in a base cylinder based on a current operation parameter The inhalation-of-air system model which 
considers that throttle opening is an orifice and calculates a throttle passage air content from a current 
throttle opening etc. is used. Based on the variation of the output of this inhalation-of-air system model, the 
variation of the restoration air content in a cylinder to intake valve close timing is predicted, this variation is 
added to the restoration air content in a base cylinder, the restoration air content in a cylinder is predicted, 
and you may make it calculate fuel oil consumption based on this restoration air content in a cylinder. Also 
in the case of a mechanical-cable-type throttle system, if it does in this way, the operation precision of the 
restoration air content in a cylinder can be raised conventionally, and the Air Fuel Ratio Control precision of 
a transient can be raised. 

[0027] Moreover, although the inhalation-of-air system model which calculates the restoration air content in 
a cylinder from a throttle opening was used in above-mentioned claims 1-14, when using the inhalation-of- 
air system model which calculates the restoration air content in a cylinder from the output (inhalation air 
content) of an intake-air-flow detection means (air flow meter), it is good like claim 1 5 for air content 
change to set the time constant of an inhalation-of-air system model as a small value also for a actual twist 
so that it may appear at an early stage. Thus, if the time constant of an inhalation-of-air system model is 
made small, the effectiveness as predicting the future restoration air content in a cylinder that a actual twist 
also has it since change of the restoration air content in a cylinder calculated with the inhalation-of-air 
system model appears at an early stage will be acquired. [ same ] Thereby, the operation precision of the 
restoration air content in a cylinder of a transient can be raised conventionally, and the Air Fuel Ratio 
Control precision of a transient can be raised. 

[0028] By the way, at claim 2 and claim 13 which were mentioned above, although the variation of the 
restoration air content in a cylinder to intake valve close timing was predicted based on the prediction 
throttle opening, since operation parameters, such as engine rotational speed, change after predicting the 
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variation of the restoration air content in a cylinder based on a prediction throttle opening before [ intake 
valve close timing ], the predictability of the restoration air content in a cylinder falls to a transient under the 
effect. 

[0029] As this cure, like claim 16, while presuming the current restoration air content in a cylinder based on 
a current throttle opening Predict a future throttle opening and the future restoration air content in a cylinder 
is predicted based on a future throttle opening. The deflection (equivalent to the prediction variation of the 
restoration air content in a cylinder) of said future restoration air content in a cylinder, and said current 
restoration air content in a cylinder It adds to the restoration air content in a base cylinder calculated based 
on the current operation parameter, the final restoration air content in a prediction cylinder is calculated, and 
you may make it calculate fuel oil consumption based on this final restoration air content in a prediction 
cylinder. If it does in this way, the prediction variation of the restoration air content in a cylinder can be 
calculated with a sufficient precision from deflection with the restoration air content in a cylinder by which 
the present was presumed to be a future restoration air content in a cylinder, and the predictability of the 
restoration air content in a cylinder can be improved. Because, the current restoration air content in a 
presumed cylinder is because information, such as the newest engine speed, is taken into consideration. (In 
claim 2 grade, it is the variation of the forecast of the restoration air content in a cylinder, and the value 
which becomes origin is calculated in the past.) 
[0030] 

[Embodiment of the Invention] «operation gestalt (1)» The operation gestalt (1) of this invention is 
hereafter explained based on drawing 1 thru/or drawing 19 . 

[003 1] First, based on drawing 1 , the outline configuration of the whole control system of the engine 1 1 
which is an internal combustion engine is explained. The upstream of the inlet pipe 12 of an engine 1 1 is 
equipped with an air cleaner 13, and the air flow meter 14 which measures an inhalation air content is 
installed in the downstream. The supply current to a heat ray is controlled so that this air flow meter 14 
keeps constant the temperature gradient of the temperature of a heat ray and the intake-air temperature 
which the heat ray (not shown) and intake-air-temperature sensing element (not shown) which are arranged 
in the flow of inhalation air are built in, and are cooled with inhalation air. According to the heat release of a 
heat ray which changes according to an intake air flow, the supply current to a heat ray changes by this, and 
the voltage signal according to this supply current is outputted as an intake-air-flow signal. 
[0032] A throttle valve 15 is formed in the downstream of this air flow meter 14, and the motors 17 (throttle 
actuator), such as a DC motor, are connected with rotation shaft 15a of this throttle valve 15. And the 
opening (throttle opening) of a throttle valve 15 is controlled by driving force of this motor 17, and this 
throttle opening is detected by the throttle opening sensor 1 8. 

[0033] In this case, a throttle opening is controlled by driving force of a motor 17, the inhalation air content 
at the time of idle operation is controlled by it also at the time of idle operation, and it carries out feedback 
control so that an engine speed may be made in agreement with target idle rotational speed. The throttle 
control at the time of this idle operation is idle rotational-speed control (ISC). In addition, idle rotational- 
speed control prepares an idle rotational-speed control bulb (ISC bulb) in the bypass path which bypasses a 
throttle valve 15, and it controls the opening of this idle rotational-speed control bulb, and you may make it 
control a bypass air content (inhalation air content at the time of idle operation) at the time of idle operation. 

[0034] On the other hand, the intake-pressure sensor 1 6 which detects an intake pressure is installed in the 
downstream of a throttle valve 15. Moreover, a fuel injection valve 20 is attached in the inlet manifold 19 
which introduces into each cylinder of an engine 1 1 the inhalation air which passed the throttle valve 15, 
and the ignition plug 21 is attached in the cylinder head of each cylinder of an engine 1 1 . The periphery of 
the signal rotor 23 attached in the crankshaft 22 of an engine 1 1 is countered, the crank angle sensor 24 is 
installed, the pulse of the engine revolution speed signal Ne outputted from this crank angle sensor 24 is 
incorporated by the electronic control unit (ECU) 25, and an engine speed is detected by the generating 
frequency of this engine revolution speed signal Ne. 

[0035] On the other hand, the amount of treading in of an accelerator pedal 26 (accelerator control input) is 
detected by the accelerator sensor 27, and the voltage signal according to this accelerator control input is 
incorporated by the electronic control unit 25 through A/D converter 28. Moreover, the output of an air flow 
meter 14, the intake-pressure sensor 16, and the various sensors of throttle opening sensor 18 grade is also 
incorporated by the electronic control unit 25 through A/D converter 28. 

[0036] It is performing the various programs for throttle control which this electronic control unit's 25 is 
constituted considering the microcomputer equipped with CPU29, ROM30, and RAM31 grade as a subject, 
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and are memorized by ROM30 by CPU29. At the time of the usual throttle control, feedback control of the 
motor 17 is carried out by PID control etc. through the motorised circuit 32 according to the opening 
command value (target throttle opening) set up based on the accelerator control input etc. A throttle opening 
is controlled by driving force of this motor 17 to an opening command value. In addition, all over the 
energization way from the motorised circuit 32 to a motor 17, the safety circuit 46 which consists of a relay 
etc. is formed, at the time of the abnormalities of an electronic throttle system, this safety circuit 46 operates 
and the energization to a motor 1 7 is intercepted. 

[0037] Furthermore, this electronic control unit 25 is performing each routine of drawing 10 memorized by 
ROM30 thru/or drawing 18 by CPU29. While performing throttle delay control mentioned later, the throttle 
opening of intake valve close timing (definite timing of the restoration air content in a cylinder) is predicted. 
The restoration air content in a cylinder is predicted based on the prediction throttle opening, fuel oil 
consumption is calculated based on the restoration air content in a prediction cylinder, the injection pulse of 
the pulse width according to the result of an operation is outputted to the actuation circuit 45, and the 
injection time (fuel oil consumption) of a fuel injection valve 20 is controlled. 

[0038] The operation approach of the fuel oil consumption by this electronic control unit 25 is explained 
based on drawing 2 thru/or drawing 9 . Drawing 2 is the block diagram showing the outline of the prediction 
approach of throttle delay control and the restoration air content in a cylinder. During engine operation, an 
accelerator control input is detected by the accelerator sensor 27, and an opening command value (target 
throttle opening) is set up by the map etc. according to an accelerator control input etc. with an opening 
command value operation means. This opening command value is predetermined time Tdly with a delay 
means. It is delayed and is outputted to the motorised circuit 32 of an electronic throttle system, time delay 
Tdly of this opening command value it is shown in drawing 3 — as — time amount Tinj from the operation 
timing (prediction timing of the restoration air content in a cylinder) of the fuel oil consumption TAU of a 
certain cylinder to the intake valve close timing of the cylinder concerned from — it is set as the time amount 
(Tdly =Tinj-Tth) which deducted the dead time Tth of an electronic throttle system. 
[0039] In this case, the dead time Tth of an electronic throttle system is the time amount Tinj from the 
operation timing (prediction timing of the restoration air content in a cylinder) of fuel oil consumption TAU 
to intake valve close timing, although it does not change even if a throttle actuation rate changes. It becomes 
so short that an engine speed becomes a high revolution. For this reason, at the time of a high revolution, it 
is the time amount Tinj from the operation timing of fuel oil consumption to intake valve close timing. 
There is a possibility of becoming shorter than the dead time Tth of an electronic throttle system. 
[0040] This point is taken into consideration and it is the time amount Tinj from the operation timing of fuel 
oil consumption TAU to intake valve close timing at this operation gestalt (1). When becoming shorter than 
the dead time Tth of an electronic throttle system, it outputs without delaying an opening command value. 
[0041] Opening command value phi total before making it delayed with a delay means on the other hand It is 
inputted into an electronic throttle model. This electronic throttle model consists of the electronic throttle 
dynamic-model section and variation operation part, as shown in drawing 4 . while this electronic throttle 
dynamic-model section simulates an electronic throttle system response delay property with a second-order 
lag element [omega2/(s2+2 zetaomegas+omega 2)] — the marginal property of the actuation rate of a throttle 
valve 1 5 — a rate limiter — simulating — opening command value phitotal before delay from — prediction 
throttle opening thetaf It calculates. Two integral elements (l-/s) of a second-order lag element are rectangle 
integrals. In addition, in order to simplify data processing, it may replace with a second-order lag element, 
and a first order lag element may be used. 

[0042] Moreover, the variation operation part of an electronic throttle model consists of a derivative element 
(d/dt) and an integral element (integral), is a derivative element (d/dt), and is the sampling time Ts of the 
output (prediction throttle opening) of the electronic throttle dynamic-model section. It asks for the 
difference of a between, and it is integrating an integral element (integral) with this difference, and 
prediction variation deltatheta of a throttle opening is calculated. Under the present circumstances, the time 
amount which integrates an integral element (integral) with difference is the time amount Tinj from the 
operation timing (prediction timing of the restoration air content in a cylinder) of fuel oil consumption TAU 
to intake valve close timing. It is the larger one of the dead time Tth of an electronic throttle system. 
Thereby, prediction variation deltatheta of the throttle opening outputted from variation operation part 
becomes the prediction variation of the throttle opening to intake valve close timing (or after dead-time Tth 
progress). 

[0043] An electronic throttle model adds prediction variation deltatheta of the throttle opening outputted 
from variation operation part to the current throttle opening theta (output of the throttle opening sensor 1 8), 
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and is prediction throttle opening thetaf of intake valve close timing (or after dead-time Tth progress). It 
asks and is this prediction throttle opening thetaf It outputs to an inhalation-of-air system model. 
[0044] As this inhalation-of-air system model is shown in drawing 5 , it consists of prediction throttle 
passage air content operation part, prediction intake-pressure operation part, and the restoration air content 
operation part in a prediction cylinder, and prediction throttle passage air content operation part considers 
that throttle opening which inhalation air passes is an orifice, and calculates the prediction throttle passage 
air content Gin from a prediction throttle opening etc. moreover, prediction intake-pressure operation part — 
the prediction throttle passage air content Gin to prediction intake pressure Pm calculating — the restoration 
air content operation part in a prediction cylinder — prediction intake pressure Pm from ~ the restoration air 
content Gcf in a prediction cylinder is calculated. Prediction throttle passage air content operation part is 
expressed with the formula of the following orifice. 
[0045] 
[Equation 3] 

Gin = it • A • Pa • f(Pm/Pa) (1) 

/R-T 

Gin : ^P^h/Vilig^CSCkg/sec] 

A : ^n^h^Mp^&fSfBaCin 2 ] 
Pa : *3?UE[Pa] 
Pm : »£Cff[Pa] 
R : ^#5£8fc 
T : RftUHJCtK] 
A = jcr 2 (l— cos 2 0 ) (2) 

Pm ^(Tfr)~^- pa »» 

f (PWPa)=y_l fr (_| r )^ --(3) 

Pm> (-TTr)'^"- pa _» 

ic : ttmt 

[0046] Here, f (Pm/Pa) is Pm/Pa, in order to simplify data processing, although you may calculate by the top 
type. It is good to compute from the table made into a parameter. The table of f (Pm/Pa) is expressed with a 
curve as shown in drawing 6 when maximum is normalized as 1 . Since this f (Pm/Pa) does not become a 
negative value physically so that clearly from the above-mentioned (3) formula and (4) types, it has been set 
as f(Pm/Pa) = 0 in the example of drawing 6 at the time of Pm/Pa>l . 

[0047] However, when it is made f(Pm/Pa) = 0 at the time of Pm/Pa >1, there is an inclination which the 
operation value (the throttle passage air content Gin, the prediction intake pressure Pm, and restoration air 
content Gcf in a prediction cylinder) of an inhalation-of-air system model vibrates, and hunting generates at 
the time of heavy load operation to which Pm/Pa is changed in the one neighborhood as shown in drawing 
8 . This cause is Pm/Pa. Since f (Pm/Pa) is guarded by 0 whenever Pm/Pa becomes one or more on an 
operation at the time of heavy load operation for the rate of change of f (Pm/Pa) to become large in the field 
of the one neighborhood, it is because change of f at the time of heavy load operation (Pm/Pa) becomes 
irregular. 
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[0048] As this cure, with this operation gestalt (1), the table of f (Pm/Pa) is set up, as shown in drawing 7 . 
That is, it is set as an f(Pm/Pa) = negative value at the time of f(Pm/Pa) =0Pm/Pa >1 at the time of f(Pm/Pa) 
= forward value Pm/Pa =1 at the time of Pm/Pa <1 . Thereby, the table of f (Pm/Pa) serves as a symmetrical 
change property which ** reverses bordering on Pm/Pa =1. 

[0049] If the table of f (Pm/Pa) of a change property like drawing 7 is used, change of f (Pm/Pa) will 
become regular at the time of heavy load operation to which Pm/Pa is changed in the one neighborhood. For 
this reason, by equalizing the operation value (the throttle passage air content Gin, the prediction intake 
pressure Pm, or restoration air content Gcf in a prediction cylinder) of an inhalation-of-air system model, as 
shown in drawing 9 , the output (restoration air content Gcf in a prediction cylinder) of the inhalation-of-air 
system model at the time of heavy load operation can be stabilized, and hunting can be prevented. 
[0050] Intake pressure Pm inputted into prediction throttle passage air content operation part Although the 
last prediction intake pressure Pm (i-1) calculated by prediction intake-pressure operation part is used, you 
may make it use the output of the intake-pressure sensor 1 6. 

[0051] Moreover, with this operation gestalt (1), although the throttle opening effective sectional area A 
used for the operation of the prediction throttle passage air content Gin may compute the throttle opening 
theta by substituting it for the aforementioned (2) formula, in order to simplify data processing, 
multiplication value mu-A of a flow coefficient mu and the throttle opening effective sectional area A is 
computed from the table which makes a prediction throttle opening a parameter. 

[0052] Next, prediction intake pressure Pm The operation approach of the restoration air content Gcf in a 
prediction cylinder is explained. If the law of conservation of mass is applied to the flow of the inhalation air 
which flows the inhalation-of-air path (henceforth a "throttle down-stream inhalation-of-air path") from a 
throttle valve 15 to the inlet port of an engine 1 1, the relation expressed with the following (5) types will be 
obtained. 

d/dt-Qm =Gin-Gcf .... (5) 

Here, it is Qm. The air content in a throttle down-stream inhalation-of-air path, and d/dt-Qm A prediction 
throttle passage air content and Gcf of the variation of the air content in a throttle down-stream inhalation- 
of-air path and Gin are the restoration air contents in a prediction cylinder. 

[0053] Moreover, if a gaseous equation of state is applied to a throttle down-stream inhalation-of-air path, 
the relation expressed to the following (6) equations will be obtained. 
Gc^eta -(Ne/2) -Vc - (Qm/VIM) .... (6) 

eta: Volumetric efficiency Ne : Engine speed Vc : Cylinder-capacity VIM: Content volume of a throttle 
down-stream inhalation-of-air path [0054] Here, volumetric efficiency eta is the engine speed Ne which is a 
parameter with an intake air flow and a correlation in order to change with intake air flows. Intake pressure 
Pm It is based and is set up on a map etc. Pm used here It is the last value Pm of a prediction intake pressure 
(i-1). 

Eta=f (Neand Pm) 

[0055] Moreover, model time constant tauIM of an inhalation-of-air system model is expressed with the 
following (7) types. 

tauIM=2andVIM/(Vc and eta-Ne ) .... (7) 

The following (8) types are drawn from the above-mentioned (5) - (7) type. 
d/dt-Qm =Gin-Qm/tauIM .... (8) 

[0056] The above-mentioned (8) formula is discretized as follows, in order to be able to carry out data 
processing of this with an electronic control unit 25, since it is continuous system. 
{Qm(i)-Qm(i-1)}/Ts =Gin(i) -Qm(i-l)/tauIM .... (9) 
Here, it is Ts. It is the sampling time. 

[0057] If this (9) type is arranged, it is the air content Qm in a throttle down-stream inhalation-of-air path. 

Operation expression is drawn as follows. 

Qm(i) ={Gin (i)-Qm(i-l)/tauIM} and Ts+Qm (i-1) [kg] 

....(10) 

[0058] moreover — if a gaseous equation of state is applied to a throttle down-stream inhalation-of-air path - 
- air content Qm in a throttle down-stream inhalation-of-air path from — prediction intake pressure Pm The 
equation to calculate is drawn as follows. 
Pm =Qm and R-T/VIM [Pa] .... (11) 

R: Gas constant T : the above-mentioned (10) formula and (11) types are used for the prediction intake- 
pressure operation part of an intake-air-temperature inhalation-of-air system model, and it is the prediction 
intake pressure Pm. It calculates. 
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[0059] The operation expression of the restoration air content Gcf in a prediction cylinder expressed with 
the following (12) types is drawn from the above-mentioned (11) formula and (6) types. 
Gcf=eta-Vc -Pm/(2 and R-T) [kg/rev] .... (12) 

The restoration air content operation part in a prediction cylinder of an inhalation-of-air system model 
calculates the temporary restoration air content Gcf in a prediction cylinder using the above-mentioned (12) 
formula. 

[0060] As shown in drawing 2 , it is inputted into a derivative element (d/dt), the difference between the 
sampling times ts is called for, and the difference integrates an integral element (integral) with the output 
(temporary restoration air content Gcf in a prediction cylinder) of an inhalation-of-air system model, the 
reset time — time amount Tinj from the operation timing (prediction timing of the restoration air content in a 
cylinder) of fuel oil consumption TAU to intake valve close timing it is . The value with which the integral 
element (integral) was integrated is prediction variation deltaGc of the restoration air content in a cylinder to 
intake valve close timing. It becomes a value considerable the bottom and is this prediction variation 
deltaGc. It adds to the restoration air content Gbase in a base cylinder calculated with the base inhalation-of- 
air system model, and the final restoration air content Gc (restoration air content in a cylinder decided to 
intake valve close timing) in a prediction cylinder is calculated. 

[0061] Next, the operation approach of the restoration air content in a base cylinder is explained. This 
restoration air content in a base cylinder is a current restoration air content in a cylinder calculated based on 
the output (intake air flow) of an air flow meter 14. Therefore, the variation of the restoration air content in a 
cylinder by change of the throttle opening from current to intake valve close timing (definite timing of the 
restoration air content in a cylinder) is not contained in this restoration air content in a base cylinder. 
Although there is an advantage that the operation precision of the restoration air content in a cylinder at the 
time of a stationary is good since the method of calculating the restoration air content in a cylinder generally 
serves as a restoration air content in an intake-air- flow = cylinder from the output of an air flow meter 14 at 
the time of a stationary Since the response delay (for example, response delay according [ the case of the air 
flow meter 14 of a heat type ] to the own heat mass of the sensor section of an air flow meter 14) of an air 
flow meter 14 exists, a transient has the fault that the responsibility of a transient is bad. 
[0062] So, with this operation gestalt (1), in order to raise the responsibility of a transient, the response 
delay of the output of an air flow meter 14 is compensated with a response delay compensator (phase-lead- 
compensation element), the output of this response delay compensator is inputted into a base inhalation-of- 
air system model, and the restoration air content Gbase in a base cylinder which is the output of this base 
inhalation-of-air system model is calculated. The transfer function of this base inhalation-of-air system 
model is expressed with the following first-order-lag type. 

Gbase=l/(l+tauIM-s), output tauIM of the restoration air-content Gdlay response delay compensator in a 
GdlayGbaseibase cylinder: Time constant [0063] Time constant tauIM of this base inhalation-of-air system 
model is expressed with a degree type. 
tauIM=2andVIM/(Vc and eta-Ne ) 

content volume Vc of the inhalation-of-air path of the VIM:throttle downstream : cylinder capacity eta: — 
volumetric efficiency Ne : an engine speed — engine speed Ne which is a parameter with an intake air flow 
and a correlation here since volumetric efficiency eta changes with intake air flows Based on an intake 
pressure P (output of the intake-pressure sensor 1 6), it is set up on a map etc. 

[0064] Prediction variation deltaGc of the restoration air content Gbase in a base cylinder calculated with 
such a base inhalation-of-air system model, and the restoration air content in a cylinder calculated from the 
prediction throttle opening etc. It integrates, the final restoration air content Gc (restoration air content in a 
cylinder decided to intake valve close timing) in a prediction cylinder is calculated, and it is this restoration 
air content Gc in a prediction cylinder. Fuel oil consumption is set up according to an engine speed etc. 
[0065] The function of each block of drawing 2 explained above is realized by each routine of drawing 10 
thru/or drawing 18 . Hereafter, the content of processing of each routine is explained to a detail. 
[0066] The main routine of [main routine] drawing 10 is performed a predetermined period after ON of an 
ignition switch. Opening command value phitotal set up according to the accelerator control input etc. when 
this routine was started, the throttle delay control routine of drawing 1 1 mentioned later was performed at 
step 1 00 and the execution condition of throttle delay control was satisfied first Predetermined time Tdly 
Throttle delay control to delay is performed. Then, the restoration air content operation routine in a 
prediction cylinder of drawing 12 which progresses to step 200 and is mentioned later is performed, and the 
restoration air content Gc (restoration air content in a cylinder decided to intake valve close timing) in a 
prediction cylinder is calculated. 
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[0067] Then, it progresses to step 300, a basic injection-quantity operation routine (not shown) is performed, 
and it is the restoration air content Gc in a prediction cylinder. Engine speed Ne It responds and is the basic 
injection quantity Tp by a map etc. It calculates. Then, it is the basic injection quantity Tp about various 
kinds of correction factors Kc, such as the fuel correction factor [ as opposed to / perform the injection- 
quantity amendment routine of drawing 18 which progresses to step 400 and is mentioned later, and / a load 
effect ] Kload (acceleration-and-deceleration correction factor), an air- fuel ratio feedback correction factor, 
and a water temperature correction factor. Multiplication is carried out and final fuel oil consumption is 
calculated. 

[0068] The throttle delay control routine of [throttle delay control-routine] drawing 1 1 is a subroutine 
performed at step 1 00 of the main routine of drawing 10 . If this routine is started, first, it responds to an 
accelerator control input (output of the accelerator sensor 27) etc. at step 101, and is opening command 
value phitotal. It sets up. under the present circumstances, opening command value phitotal Demand 
opening phipedal according to an accelerator control input Demand opening phiisc by idle rotational-speed 
control (ISC) etc. — it integrates and asks for various demand openings. 

phitotal =phipedal+phiisc [0069] Processing of this step 101 plays a role of an opening command value 
operation means as used in the field of a claim. 

[0070] Then, it progresses to step 1 02 and judges whether throttle delay control prohibition conditions are 
satisfied, it is that the time of that it is predetermined within a time [ immediately after the time of ** start 
up or start up ] as throttle delay control prohibition conditions here, for example and ** idle operation or an 
accelerator control input is small, that ** automatic transmission is in a neutral condition, etc., if there are 
conditions which correspond at least one of these conditions, throttle delay control prohibition conditions 
will be satisfied, and when other, throttle delay control prohibition conditions become abortive. 
[0071] The time of start up, and immediately after start up, since the engine revolution is unstable, when 
throttle delay control is performed essentially, there is a possibility that engine revolution fluctuation may 
become still larger. Moreover, at the time of idle operation, in order for idle rotational-speed control (ISC) to 
operate and to carry out feedback control of the idle rotational speed, when throttle delay control is 
performed, there is a possibility that idle rotational-speed control may interfere with throttle delay control, 
and an idle revolution may become instability. Moreover, since an operator may perform racing (an engine 
carrying out a fore blow) when an automatic transmission is in a neutral condition, when it is in a neutral 
condition and throttle delay control is performed, the start of an engine speed is overdue at the time of racing 
activation, and there is a possibility that accelerator responsibility and acceleration nature may impress an 
operator with it being bad. 

[0072] So, with this operation gestalt (1), the adverse effect by throttle delay control is lost by forbidding 
throttle delay control, when it is the operational status in which the adverse effect by throttle delay control 
appears (at the time [ At the at the time of idle operation immediately after start up ] of a neutral). [ At the 
time of start up ] 

[0073] If judged with throttle delay control prohibition conditions being satisfied at step 102, throttle delay 
control will be forbidden, and it will progress to step 103, and will output to the motorised circuit 32, 
without delaying current (newest) opening command value phitotal (i). 

[0074] On the other hand, at step 102, if judged with throttle delay control prohibition conditions being 
abortive, throttle delay control will be carried out as follows by processing after step 104. First, it is opening 
command value phitotal at step 1 04. Time delay Tdly It determines, under the present circumstances, time 
delay Tdly it is shown in drawing 3 — as — time amount Tinj from the operation timing (prediction timing of 
the restoration air content in a cylinder) of fuel oil consumption TAU to intake valve close timing from » it 
is set as the time amount (Tdly =Tinj-Tth) which deducted the dead time Tth of an electronic throttle 
system. However, time amount Tinj from the operation timing of fuel oil consumption TAU to intake valve 
close timing When becoming shorter than a dead time Tth (in the case of Tinj-Tth<0), it is a time delay 
Tdly. It is referred to as 0. 

[0075] Then, it progresses to step 105 and is a time delay Tdly. The inner number Cdly of samplings It 
calculates by the degree type. 

Cdly =Tdly / Ts — here — Ts It is the sampling time. 

[0076] Then, it progresses to step 106 and is a time delay Tdly from current. The inner number Cdly of 
samplings Opening command value phitotal (i-Cdly) calculated before is outputted to the motorised circuit 
32. thereby — opening command value phitotal output timing — time delay Tdly only — you make it delayed 
Processing of these steps 102-106 plays a role of a delay means as used in the field of a claim. 
[0077] The restoration air content operation routine in a prediction cylinder of [restoration air content 
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operation routine in prediction cylinder] drawin g 12 R> 2 is a subroutine performed at step 200 of the main 
routine of drawing 11 , and plays a role of a restoration air content operation means in a prediction cylinder 
as used in the field of a claim. 

[0078] If this routine is started, first, at step 201, the prediction intake-pressure operation routine of drawing 

13 mentioned later will be performed, and the prediction intake pressure Pm (intake pressure of intake valve 
close timing) will be calculated. Then, it progresses to step 202 and is the prediction intake pressure Pm. It 
uses and is the restoration air content Gcf in a prediction cylinder (i) by the degree type. It calculates. 

Gcf (i) = eta-Vc -Pm/(2 and R-T) [kg/rev] 

eta: Volumetric efficiency Vc : Cylinder-capacity Rrgas constant T: Intake-air temperature [0079] Then, it 
progresses to step 203 and judges whether throttle delay control prohibition conditions are satisfied by the 
same approach as step 102 of drawing 1 1 . It progresses to step 204 from step 203, without carrying out 
throttle delay control, when throttle delay control prohibition conditions are satisfied, and is prediction 
variation deltaGc of the restoration air content in a cylinder from the operation timing of fuel oil 
consumption to intake valve close timing. It is set as 0. 

[0080] On the other hand, when throttle delay control prohibition conditions are abortive, throttle delay 
control is carried out by the throttle delay control routine of drawing 1 1 , and it progresses to step 205 from 
step 203, and is the prediction time amount Tinj. The inner number Cp of samplings It calculates by the 
degree type. 

Cp =Tinj / Ts — here ~ prediction time amount Tinj the time amount from the operation timing of fuel oil 
consumption to intake valve close timing - it is - Ts It is the sampling time. 

[0081] Then, it progresses to step 206 and is prediction variation deltaGc of the restoration air content Gcf 
in a cylinder from the operation timing of fuel oil consumption to intake valve close timing. It calculates by 
the degree type. 

deltaGc = - Gcf (i) -Gcf (i-Cp) - here Gcf (i) this restoration air content in a prediction cylinder (that is, 
restoration air content in a prediction cylinder of the operation timing of intake valve close timing) - it is - 
Gcf (i-Cp) - current - prediction time amount Tinj The inner number Cp of samplings It is the restoration 
air content in a cylinder (that is, restoration air content in a cylinder of the operation timing of fuel oil 
consumption) calculated before. 

[0082] Prediction variation deltaGc It progresses to step 207 after an operation, the restoration air content 
operation routine in a base cylinder (not shown) is performed, and the restoration air content Gbase in a base 
cylinder is calculated. Under the present circumstances, the response delay of the output of an air flow meter 

14 is compensated with a response delay compensator (phase-lead-compensation element), and the 
restoration air content Gbase in a base cylinder is calculated with the following transfer function using the 
output Gdlay of this response delay compensator. 

Gbase=l/and (1+tauIM-s), and Gdlay - here, tauIM is a time constant. Although the operation expression of 
the restoration air content in a base cylinder is expressed with the continuous system in order that a top type 
may simplify explanation, an electronic control unit 25 calculates the restoration air content Gbase in a base 
cylinder using the formula discretized in the top type. 

[0083] Then, prediction variation deltaGc for which progressed to step 208 and this restoration air content 
Gbase in a base cylinder was asked at the above-mentioned step 206 It adds and is the final restoration air 
content Gc in a prediction cylinder. It asks. 

Gc = Gbase+delta Gc [0084] The prediction intake-pressure operation routine of [prediction intake-pressure 
operation routine] drawing 13 is a subroutine performed at step 201 of the restoration air content operation 
routine in a prediction cylinder of drawing 12 . If this routine is started, first, at step 21 1, the prediction 
throttle passage air content operation routine of drawing 14 mentioned later will be performed, and the 
prediction throttle passage air content Gin will be calculated. Then, the inhalation-of-air system model time 
constant operation routine of drawing 16 R> 6 which progresses to step 212 and is mentioned later is 
performed, and model time constant tauIM of an inhalation-of-air system model is calculated. Then, it 
progresses to step 213 and is the air content Qm in a throttle down-stream inhalation-of-air path. It 
calculates by the degree type. 

Qm(i) ={Gin (i)-Qm(i-l)/tauIM} and Ts+Qm (i-1) - here - Qm (i) - the air content in this throttle down- 
stream inhalation-of-air path, the air content in the throttle down-stream inhalation-of-air path of last time 
[ Qm / (i-1) ], and Ts It is the sampling time. 

[0085] then, the step 214 - progressing ~ air content Qm in a throttle down-stream inhalation-of-air path 
from — prediction intake pressure Pm It calculates by the degree type. 

Pm =Qm and R-T/VIM - here, R is [ an intake-air temperature and VIM of a gas constant and T ] the 
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content volume of a throttle down-stream inhalation-of-air path. 

[0086] Then, it is the prediction intake pressure Pm by progressing to step 215 and computing the average of 
this prediction intake pressure Pm (i) and the last prediction intake pressure Pm (i-1). It equalizes. 
Pm (i) = {Pm(i)+Pm (i-l)}/2 [0087] The prediction throttle passage air content operation routine of 
[prediction throttle passage air content operation routine] drawing 14 is a subroutine performed at step 211 
of the prediction intake-pressure operation routine of drawing 13 . If this routine is started, the prediction 
throttle opening operation routine of drawing 15 mentioned later is first performed at step 221 , and it is 
prediction throttle opening thetaf of intake valve close timing. It calculates. Then, it progresses to step 222 
and atmospheric pressure Pa, an intake-air temperature T, and the last prediction intake pressure Pm (i-1) are 
read. 

[0088] Then, it progresses to step 223 and the prediction throttle passage air content Gin is calculated by the 
degree type. 
[0089] 
[Equation 4] 

Gin = » • A • Pa • f (Pm/Pa) 
/R-T 

Gin : ^P^h/loliSS^ftCkg/sec] 

* : m&&& 

Pa : *3cCJBE[Pa] 

Pm : »«ff[Pa] 
R : 

T : VMSMUtOQ 

f(Pm/Pa) : KtePmt^^Pa<Dlfc-e^*5^Sffi 

[0090] Under the present circumstances, mu-A is prediction throttle opening thetaf. Computing from the 
table made into a parameter, f (Pm/Pa) is Pm/Pa. It computes from the table of drawing 7 made into a 
parameter. Intake pressure Pm The last prediction intake pressure Pm (i-1) is used, and, as for atmospheric 
pressure Pa and an intake-air temperature T, the detection value of a sensor is used, respectively. In addition, 
atmospheric pressure Pa A standard atmosphere (fixed value) may be used. 

[0091] The prediction throttle opening operation routine of [prediction throttle opening operation routine] 
drawing 1 5 R> 5 is a subroutine performed at step 221 of the prediction throttle passage air content 
operation routine of drawing 14 , and plays a role of a prediction throttle opening operation as used in the 
field of a claim. 

[0092] If this routine is started, first, it responds to an accelerator control input etc. at step 231, and is 
opening command value phitotal. It sets up. under the present circumstances, opening command value 
phitotal Demand opening phipedal according to an accelerator control input Demand opening phiisc by idle 
rotational-speed control (ISC) etc. — it integrates and asks for various demand openings, 
phitotal =phipedaH-phiisc [0093] Then, it is opening command value phitotal before delay by the electronic 
throttle dynamic-model section and variation operation part of the electronic throttle model which 
progresses to step 232, progresses to step 233 after reading the present throttle opening theta detected by the 
throttle opening sensor 18, and is shown in drawing 4 . It uses and prediction variation deltatheta of a 
throttle opening is calculated. This prediction variation deltatheta is the time amount Tinj from the operation 
timing (prediction timing of the restoration air content in a cylinder) of fuel oil consumption TAU to intake 
valve close timing. It is the prediction variation of a throttle opening. However, time amount Tinj to intake 
valve close timing When shorter than the dead time Tth of an electronic throttle system, prediction variation 
deltatheta of the throttle opening within this dead time Tth is called for. 

[0094] Then, it progresses to step 234, prediction variation deltatheta is added to the current throttle opening 
theta, and it is prediction throttle opening thetaf. It asks. 

thetaf = theta+delta theta — this prediction throttle opening thetaf It is the prediction throttle opening of 
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intake valve close timing (or after dead-time Tth progress). 

[0095] The inhalation-of-air system model time constant operation routine of [inhalation-of-air system 
model time constant operation routine] drawing 16 R> 6 is a subroutine performed at step 212 of the 
prediction intake-pressure operation routine of drawing 13 . If this routine is started, first, at step 241, the 
volumetric-efficiency operation routine of drawing 1 7 mentioned later will be performed, and volumetric 
efficiency eta will be calculated. Then, it progresses to step 242 and model time constant tauIM is calculated 
by the degree type. 

tauIM=2andVIM/(eta-Ne [ VC and ]/60) 

Here, VIM is the content volume (fixed value) of a throttle down-stream inhalation-of-air path, and VC. 
Cylinder capacity (fixed value) and Ne It is an engine speed (rpm). 

[0096] The volumetric-efficiency operation routine of [volumetric-efficiency operation routine] drawing 17 
is a subroutine performed at step 241 of the inhalation-of-air system model time constant operation routine 
of drawing 16 . Starting of this routine first reads the last intake pressure Pm (i-1) and the last atmospheric 
pressure Pa, an intake-air temperature T, an engine speed Ne, valve timing WT, and the cooling water 
temperature THW at step 151. Then, basic volumetric-efficiency etar progress to step 152, search the 
volumetric-efficiency map which makes a parameter Pm/Pa, an engine speed Ne, and valve timing WT, 
and corresponding to the current engine operation condition It calculates and is this basic volumetric- 
efficiency etar. It amends with the correction value according to the cooling water temperature THW, and 
volumetric efficiency eta is searched for. 

[0097] The injection-quantity amendment routine of [injection-quantity amendment routine] drawing 18 is a 
subroutine performed at step 400 of the main routine of drawing 10 , and plays a role of a fuel-oil- 
consumption operation means as used in the field of [ routine / (not shown) / basic injection-quantity 
operation ] a claim. 

[0098] a ****** [ whether it is a load effect (fluctuation of the restoration air content in a cylinder) by 
accelerator actuation at step 401, and that for example, an accelerator control input is beyond the set point 
first when this routine is started ] ~ or it judges by whether the variation of an accelerator control input is 
beyond the set point. It is the small value Kl about the fuel correction factor [ as opposed to / if judged with 
the load effect by accelerator actuation / progress to step 402 and / a load effect (fluctuation of the 
restoration air content in a cylinder) ] Kload. It sets up. By the operation approach of the restoration air 
content in a cylinder of this operation gestalt (1), since this reason can predict the load effect (fluctuation of 
the restoration air content in a cylinder) by accelerator actuation with a sufficient precision, it is because 
amendment to fuel oil consumption can be lessened. 

[0099] It is the large value K2 about the fuel correction factor [ as opposed to / when judged with it not 
being a load effect by accelerator actuation (in for example, the case of the load effect by the time of shifting 
an automatic transmission from a neutral range to a drive range, power steering, a brake, an air-conditioner, 
etc.) / progress to step 403 and / a load effect on the other hand ] Kload. It sets up. Since this reason cannot 
predict the load effect by factors other than accelerator actuation from an accelerator control input, it is 
because it is desirable to make [ many ] amendment to fuel oil consumption to the load effect by factors 
other than accelerator actuation. 

[0100] After determining the fuel correction factor [ as opposed to a load effect as mentioned above ] Kload 
by step 402 or 403, it progresses to step 404 and various kinds of fuel correction factors Kc to factors other 
than a load effect (for example, an air- fuel ratio feedback correction factor, a water temperature correction 
factor, a study correction factor, etc.) are calculated, and they are the basic injection quantity Tp and the fuel 
correction factors Kload and Kc at the following step 405. And invalid injection time Tv It uses and the final 
fuel oil consumption (injection pulse width) TAU is calculated by the degree type. 

TAU=Tp xKloadxKc+Tv[0101] An example of the behavior of the prediction throttle opening calculated by 
each routine explained above and the restoration air content in a prediction cylinder is shown in the timing 
diagram of drawing 19 . under engine operation — an accelerator control input etc. — responding — opening 
command value phitotal setting up — this opening command value phitotal output timing — time delay Tdly 
only — you make it delayed under the present circumstances, time delay Tdly it is shown in drawing 3 — as - 
- time amount Tinj from the operation timing (prediction timing of the restoration air content in a cylinder) 
of fuel oil consumption TAU to intake valve close timing from — it is set as the time amount (Tdly =Tinj- 
Tth) which deducted the dead time Tth of an electronic throttle system. However, time amount Tinj from the 
operation timing of fuel oil consumption TAU to intake valve close timing When becoming shorter than the 
dead time Tth of an electronic throttle system (in the case of Tinj-Tth<0), it is a time delay Tdly. It is 
referred to as 0. 
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[0102] Opening command value phi total before delay It is based, prediction variation deltatheta of a throttle 
opening is calculated with the electronic throttle model of drawing 4 R> 4, this prediction variation 
deltatheta is added to the present throttle opening theta (output of the throttle opening sensor 18), and it is 
prediction throttle opening thetaf of intake valve close timing (or after dead-time Tth progress). It asks. And 
this prediction throttle opening thetaf It uses, the temporary restoration air content Gcf in a prediction 
cylinder is calculated with the inhalation-of-air system model of drawing 5 , differential / integral processing 
of this is carried out, and it is prediction variation deltaGc of the restoration air content in a cylinder to 
intake valve close timing. It calculates. This prediction variation deltaGc It adds to the restoration air content 
Gbase in a base cylinder calculated with the base inhalation-of-air system model, and the final restoration air 
content Gc (restoration air content in a cylinder decided to intake valve close timing) in a prediction cylinder 
is calculated. Thereby, it is the restoration air content Gc in a cylinder. It can become possible to predict 
with a sufficient precision, and the Air Fuel Ratio Control precision of a transient can be raised. 
[0103] «operation gestalt (2)» With the above-mentioned operation gestalt (1), although the delay means 
(refer to drawing 2 ) was made to perform throttle delay control, in the operation gestalt (2) of this invention 
shown in drawing 20 , a delay means is excluded, and throttle delay control is not performed, but a throttle 
opening is predicted using the dead time Tth of an electronic throttle system. 

[0104] With this operation gestalt (2), it outputs to the motorised circuit 32, without delaying the opening 
command value set up based on the accelerator control input etc. with the opening command value operation 
means. By the same approach as said operation gestalt (1), and with an electronic throttle model Based on an 
opening command value and the present throttle opening (output of the throttle opening sensor 1 8), the 
throttle opening of intake valve close timing (or after dead-time Tth progress) is predicted. With an 
inhalation-of-air system model (configuration of drawing 5 ) The temporary restoration air content in a 
prediction cylinder is calculated from this prediction throttle opening, differential / integral processing of 
this is carried out, and the prediction variation of the restoration air content in a cylinder to intake valve 
close timing (or after dead-time Tth progress) is calculated. And it adds to the restoration air content in a 
base cylinder which calculated this prediction variation with the base inhalation-of-air system model, and 
the final restoration air content in a prediction cylinder is calculated. 

[0105] Also in this operation gestalt (2) explained above, it can become possible to predict a throttle opening 
using the dead time Tth of an electronic throttle system, and to predict the restoration air content in a 
cylinder with a sufficient precision from the prediction throttle opening, and the Air Fuel Ratio Control 
precision of a transient can be raised. 

[0106] «operation gestalt (3)» Although the above-mentioned operation gestalt (1) and (2) apply this 
invention to an engine with an electronic throttle system, the operation gestalt (3) shown in drawing 21 
applies this invention to the engine with a mechanical-cable-type throttle system which interlocks a throttle 
opening with accelerator actuation mechanically. 

[0107] With this operation gestalt (3), since an accelerator control input and a real throttle opening interlock 
mechanically and there is no response delay of a throttle, the opening command value operation means, the 
delay means, and the electronic throttle model are not prepared. Although the prediction throttle opening 
was inputted into the inhalation-of-air system model, with this operation gestalt (3), a current throttle 
opening (output of the throttle opening sensor 18) is inputted into an inhalation-of-air system model said 
operation gestalt (1) and (2). The configuration of this inhalation-of-air system model is substantially [ as 
said operation gestalt (1) ] the same, calculates the temporary restoration air content in a prediction cylinder 
from a current throttle opening, carries out differential / integral processing of this, and calculates the 
prediction variation of the restoration air content in a cylinder to intake valve close timing (or after 
predetermined period progress). And it adds to the restoration air content in a base cylinder which calculated 
this prediction variation with the base inhalation-of-air system model, and the final restoration air content in 
a prediction cylinder is calculated. 

[0108] Also in the case of a mechanical-cable-type throttle system, if it does in this way, the operation 
precision of the restoration air content in a cylinder can be raised conventionally, and the Air Fuel Ratio 
Control precision of a transient can be raised. 

[0109] «operation gestalt (4)» in said operation gestalt (1) and (2) As shown in drawing 2 and drawing 
20 , based on an opening command value and the present throttle opening (output of the throttle opening 
sensor 1 8), the throttle opening of intake valve close timing (or after dead-time Tth progress) is predicted 
with an electronic throttle model. With an inhalation-of-air system model Calculate the temporary 
restoration air content in a prediction cylinder from this prediction throttle opening, and differential / 
integral processing of this is carried out. After calculating the prediction variation of the restoration air 
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content in a cylinder to intake valve close timing (or after dead-time Tth progress), he adds to the restoration 
air content in a base cylinder which calculated this prediction variation with the base inhalation-of-air 
system model, and is trying to calculate the final restoration air content in a prediction cylinder. 
[01 10] on the other hand, with the operation gestalt (4) of this invention shown in drawing 22 thru/or 
drawing 25 With an electronic throttle model, the throttle opening of intake valve close timing (or after 
dead-time Tth progress) is predicted based on an opening command value and a current throttle opening 
(output of the throttle opening sensor 18). While calculating the future restoration air content in a cylinder 
(temporary restoration air content in a prediction cylinder) from this prediction throttle opening, with an 
inhalation-of-air system model with an inhalation-of-air system model Based on a current throttle opening 
(output of the throttle opening sensor 1 8), the current restoration air content in a cylinder is calculated. The 
deflection (equivalent to the prediction variation of the restoration air content in a cylinder) of the future 
restoration air content in a cylinder, and the current restoration air content in a cylinder It adds to the 
restoration air content in a base cylinder calculated with the base inhalation-of-air system model, and he 
calculates the final restoration air content in a prediction cylinder, and is trying to calculate fuel oil 
consumption based on this final restoration air content in a prediction cylinder. 
[0111] Except for the point of performing the same main routine as drawing 10 , and performing the 
restoration air content operation routine in a prediction cylinder of drawing 23 at step 200 also by this 
operation gestalt (4), it is the same as said operation gestalt (1). 

[0112] By the restoration air content operation routine in a prediction cylinder of drawing 23 , the present 
restoration air content presumption routine in a cylinder of drawing 24 mentioned later is performed by step 
500, and it is based on the present throttle opening theta (output of the throttle opening sensor 1 8) with an 
inhalation-of-air system model, and is the present restoration air content Gest in a cylinder. It calculates. 
[0113] then, the future restoration air content operation routine in a cylinder of drawing 25 which progresses 
to step 600 and is mentioned later — performing — an electronic throttle model — an opening command value 
and the present throttle opening theta — being based — throttle opening thetaf of intake valve close timing 
(or after dead-time Tth progress) predicting - an inhalation-of-air system model — this prediction throttle 
opening thetaf from — the future restoration air content Gcf (temporary restoration air content in a prediction 
cylinder) in a cylinder is calculated. 

[0114] Then, restoration air content Gest in a cylinder progress, the future restoration air content Gcf in a 
cylinder, and current to step 700 after progressing to step 600 and calculating the restoration air content 
Gbase in a base cylinder by the same approach as said operation gestalt (1) Deflection (equivalent to the 
prediction variation of the restoration air content in a cylinder) is added to the restoration air content Gbase 
in a base cylinder, and it is the final restoration air content Gc in a prediction cylinder. It asks. 
Gc =Gbase+(Gcf-Gest ) 

[0115] In the present restoration air content presumption routine in a cylinder of drawing 24 , the present 
throttle opening theta is first read at step 501, and they are atmospheric pressure Pa, an intake-air 
temperature T, and an intake pressure Pm at the following step 502. It reads. Under the present 
circumstances, intake pressure Pm The detection value of the intake-pressure sensor 1 6 may be used, or the 
last value of the prediction intake pressure calculated at step 601 of drawing 25 mentioned later may be 
used. 

[0116] Then, it progresses to step 503 and the present throttle passage air content Gin is calculated by the 
same approach as the routine of drawing 14 explained with said operation gestalt (1). Then, it progresses to 
step 504 and model time constant tauIM of an inhalation-of-air system model is calculated by the same 
approach as the routine of drawing 16 explained with said operation gestalt (1). 

[0117] Then, it is the air content Qm in a throttle down-stream inhalation-of-air path with the same approach 
as step 213 of drawing 13 which progressed to step 505 and was explained with said operation gestalt (1). It 
calculates by the degree type. 

Qm(i) ={Gin (i)-Qm(i-l)/tauIM} and Ts+Qm (i-1) - here ~ Qm (i) - the air content in this throttle down- 
stream inhalation-of-air path, the air content in the throttle down-stream inhalation-of-air path of last time 
[ Qm / (i-1) ], and Ts It is the sampling time. 

[0118] then, the step 506 — progressing — air content Qm in a throttle down-stream inhalation-of-air path 
from ~ current intake pressure Pm It calculates by the degree type. 

Pm =Qm and R-T/VIM - here, R is [ an intake-air temperature and VIM of a gas constant and T ] the 
content volume of a throttle down-stream inhalation-of-air path. 

[0119] Then, it is an intake pressure Pm by progressing to step 507 and computing the average of this intake 
pressure Pm (i) and the last intake pressure Pm (i-1). It equalizes. 
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Pm (i) = {Pm(i)+Pm (i-l)}/2 [0120] Then, it progresses to step 508 and is this intake pressure Pm. It uses 
and is the current restoration air content Gest in a cylinder by the degree type. It calculates. 
Gest = eta-Vc -Pm/(2 and R-T) 

eta: Volumetric efficiency Vc : Cylinder-capacity R:gas constant T: Intake-air temperature [0121] On the 
other hand, in the future restoration air content operation routine in a cylinder of drawing 25 , the prediction 
intake pressure Pm (intake pressure of intake valve close timing) is first calculated at step 601 by the same 
processing as the prediction intake-pressure operation routine of drawing 13 explained with said operation 
gestalt (1). Then, it progresses to step 602 and is the prediction intake pressure Pm. It uses and the future 
restoration air content Gcf (restoration air content in a cylinder of intake valve close timing) in a cylinder is 
calculated by the degree type. 
Gcf^eta-Vc -Pm/(2 and R-T) 

[0122] With this operation gestalt (4) explained above, while presuming the current restoration air content in 
a cylinder based on a current throttle opening Since a future throttle opening is predicted, the future 
restoration air content in a cylinder is predicted and the prediction variation of the restoration air content in a 
cylinder was calculated from the deflection of said future restoration air content in a cylinder, and said 
current restoration air content in a cylinder From said operation gestalt (1) and (2), the prediction variation 
of the restoration air content in a cylinder can be calculated with a sufficient precision, and the predictability 
of the restoration air content in a cylinder can be improved. 

[0123] «operation gestalt (5)» Although the inhalation-of-air system model which calculates the 
restoration air content in a cylinder from a throttle opening was used in above-mentioned operation gestalt 
(1) - (4), with the operation gestalt (5) of this invention, air content change sets the time constant of this 
inhalation-of-air system model as a small value also for a actual twist using the inhalation-of-air system 
model which calculates the restoration air content in a cylinder from the output (intake air flow) of an air 
flow meter 14 (intake-air-flow detection means) so that it may appear at an early stage. 
[0124] With this operation gestalt (5), the time constant of an inhalation-of-air system model is set as a 
small value instead of establishing a means to predict the restoration air content in a cylinder. If it does in 
this way, the effectiveness as predicting the future restoration air content in a cylinder that a actual twist also 
has it since change of the restoration air content in a cylinder calculated with the inhalation-of-air system 
model appears at an early stage will be acquired. [ same ] Thereby, the operation precision of the restoration 
air content in a cylinder of a transient can be raised conventionally, and the Air Fuel Ratio Control precision 
of a transient can be raised. 
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* NOTICES * 

JPO and NCI PI are not responsible for any 
damages caused by the use of this translation. 

1 .This document has been translated by computer. So the translation may not reflect the original precisely. 
2.**** shows the word which can not be translated. 
3. In the drawings, any words are not translated. 
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